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ABSTRACT 



The physics of a twisted magnetic flux tube rising in a stratified medium is studied 
using a numerical MHD code. The problem considered is fully compressible (no Boussi- 
nesq approximation), includes ohmic resistivity, and is two dimensional, i.e., there is no 
variation of the variables in the direction of the tube axis. We study a high plasma (3 
case with small ratio of radius to external pressure scaleheight. The results obtained can 
therefore be of relevance to understand the transport of magnetic flux across the solar 
convection zone. 

'■ We confirm that a sufficient twist of the field lines around the tube axis can suppress 

the conversion of the tube into two vortex rolls. For a tube with relative density deficit 
of order 1/(3 (the classical Parker buoyancy) and radius smaller than a pressure scale- 
XJ : height (i?2 < i//), the minimum amount of twist necessary corresponds to an average 

pitch angle of order sin~^[(i?/i7p)^/^]. The evolution of a tube with this degree of twist is 
studied in detail, including the initial transient phase, the internal torsional oscillations 
and the asymptotic, quasi-stationary phase. During the initial phase, the outermost, 
weakly magnetized layers of the tube are torn off its main body and endowed with 
vorticity. They yield a trailing magnetized wake with two vortex rolls. Which fraction 
of the total magnetic flux is brought to the wake is a function of the initial degree of 
twist. In the weakly twisted case, most of the initial tube is turned into vortex rolls. 
At the opposite end (strong initial twist), the tube rises with only a small deformation 
and no substantial loss of magnetic flux. The formation of the wake and the loss of flux 
from the main body of the tube are basically complete after the initial transient phase. 

A sharp interface between the tube interior and the external flows is formed at the 
tube front and sides; it has the characteristic features of a magnetic boundary layer. Its 
structure is determined as an equilibrium between ohmic diffusion and fleld advection 
through the external flows. It is the site of vorticity generation via the magnetic fleld 
during the whole tube evolution. 
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From the hydrodynamical point of view, this problem constitutes an intermediate 
case between the rise of air bubbles in water and the motion of a rigid cylinder in an 
external medium. As in the first one, the tube is deformable and the outcome of the 
experiment (shape of rising object and wake) depends on the value of the Weber number. 
Several structural features obtained in the present simulation arc also observed in rising 
air bubbles, like a central tail, and a skirt enveloping the wake. As in rigid cylinders, the 
boundary layer satisfies a no-slip condition (provided for in the tube by the magnetic 
field), and secondary rolls are formed at the lateral edges of the moving object. 

Subject headings: Magnetohydrodynamics: MHD, Hydrodynamics, Sun: Magnetic 
Fields, Magnetic Fields, Sun: Activity, Stars: Magnetic Fields 
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1. Introduction 

The rise of magnetic flux from the deep levels of 
the solar eonvection zone to the photosphere is a 
complex phenomenon involving many different mag- 
netie and hydrodynamical processes. Particular at- 
tention has been devoted in the past 20 years to 
the time-evolution of a single buoyant magnetic flux 
tube considered as a one-dimensional object (Moreno- 
Insertis 1986; Choudhuri 1989; D'Silva & Choudhuri 
1993; Fan, Fisher & DeLuca 1993; Fan, Fisher & Mc- 
Clymont 1994; Caligari et al. 1995. Further refer- 
ences and a recent review can be found in the paper 
by Moreno-Insertis 1997a). These calculations incor- 
porate several aspects of the basic physics of the rise 
of the magnetic tubes (buoyancy, magnetic and rota- 
tional forces, external stratification, etc). They have 
been successful in predicting morphological and kine- 
matic features of the resulting active regions which 
are observed at the surface of the Sun. 

Yet, the assumption of one-dimensionality of the 
magnetic region is certainly a drastic simplification. 
From laboratory and numerical experiments in dif- 
ferent contexts we know of the complicated hydro- 
dynamical and magnetic structure within and out- 
side a tubular object which is moving with respect to 
the surrounding fluid. The one-dimensional numer- 
ical models mentioned above, in particular, do not 
contain two ingredients which turn out to be funda- 
mental in studying the rise of buoyant magnetized 
plasma regions, to wit, the vorticity of the velocity 
field and the twist of the field lines around the main 
axis of the tube. 

Vorticity and transverse field components may be 
crucial for the formation of the tubes in the first 
place (Cattaneo & Hughes 1988; Cattaneo, Chiueh, 
& Hughes 1990; Matthews, Hughes & Proctor 1995). 
They also play a central role in the time evolution 
of the rising magnetic region. A clear warning in 
this sense came from the work of Schiissler (1979), 
who showed how the cross section of a straight, buoy- 
ant magnetic tube initially with the same tempera- 
ture of its surroundings develops an umbrella shape 
(two side lobes connected on their upper side by an 
arch). The side lobes rotate in opposite directions 
around a horizontal axis, each thus constituting a 
vortex tube; they finally detach from each other and 
from the arch above them. The whole process occurs 
at the beginning of the rise, namely before the tube 
has risen across a height equivalent to a few times its 



own diameter. The physics involved has been consid- 
ered in detail by Longcope, Fisher & Arendt (1996). 
These authors have studied the Boussinesq problem, 
including untwisted and very weakly twisted magnetic 
tubes. They clearly show how the two rotating side 
lobes, when detached from the rest, are subjected to a 
downward-pointing lift force, as a result of their flow 
being non-circulation free. The lift ends up cancelling 
the buoyancy force, this being the reason for their 
horizontal asymptotic motion. If this were a univer- 
sal mechanism operating on all rising magnetic flux 
tubes, then magnetic buoyancy should no longer be 
considered an efficient mechanism to bring magnetic 
flux to the photosphere. 

In the present paper we consider in detail the more 
general case of a buoyant magnetic flux tube with 
an arbitrary initial twist of the field lines (but still 
horizontal and with uniform values of all variables 
along the direction of the axis) . The transverse mag- 
netic field (i.e., the component of the field vector nor- 
mal to the tube axis) imparts a certain rigidity to 
the tube cross section. If strong enough, it can pre- 
vent the conversion of the tube into a vortex tube 
pair. The minimum amount of twist necessary for 
that corresponds to an average pitch angle of order 
sm-^{[R/Hp (3/2 |Ap/p|]V2}, with R the tube ra- 
dius. Hp the external pressure scaleheight, and Ap 
the density difference between tube and surround- 
ings (see §3.1). This approximate criterion is indeed 
fulfilled by magnetic tubes with the classical Parker 
magnetic buoyancy, as shown in a preliminary pre- 
sentation of this paper (Moreno-Insertis & Emonet 
1996). In the present article, we explore in more de- 
tail the physics involved in that process, discussing 
a number of (M)HD processes occurring inside the 
tube, in the boundary layer at its periphery and in 
the trailing wake. We also show how the results of 
Schiissler (1979) and Longcope et al.(1996) can be 
seen as the limiting case in which the trailing wake 
in fact engulfs most of the original rising tube. The 
pitch angle just mentioned thus signals the borderline 
between the weak and strong twist regimes: a buoy- 
ant tube with an initial twist above that level rises 
without being strongly deformed and is followed by 
a wake containing only a small fraction of the initial 
total magnetic flux. 

Additionally to the work of Moreno-Insertis & 
Emonet (1996), there is another paper in the recent 
literature dealing with a related subject (Cargill et 
al. 1996). The authors have studied the interaction 
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of a twisted tube with a magnetized medium in the 
absence of gravity when the tube is subjected to an 
ad- hoc, spatially uniform acceleration. Special em- 
phasis was put on the reconnection of the ambient 
magnetic field with the tube's own one. Buoyancy, 
stratification, or diff'erent degrees of twist were not 
studied in that paper. Two further papers, submit- 
ted simultaneously with the present one, deal with the 
rise of buoyant twisted magnetic tubes (Fan, Zweibel 
& Lantz 1997; Hughes, Falle & Joarder, 1998). The 
first authors, in particular, study the interaction be- 
tween tubes rising in pairs. The results of both pa- 
pers concerning the rise of single tubes are in general 
agreement with those of Moreno-Insertis & Emonet 
(1996). 

The layout of the paper is as follows. After a brief 
presentation of the equations and the numerical pro- 
cedure in §2, the basic features of the physical prob- 
lem are considered (§3). This includes the main pa- 
rameters and a discussion of the amount of twist nec- 
essary to prevent the deformation of the tube and its 
conversion into vortex rolls. In §4, the simulation of 
the rise of a tube with that amount of twist is pre- 
sented. In particular, the initial acceleration phase, 
the internal torsional oscillations and the later asymp- 
totic phase are discussed. Section §5 deals with the 
structure of the magnetic boundary layer around the 
tube, and §6 examines the trailing wake. Finally, the 
transition between the twisted and untwisted case is 
explained in §7. A general discussion follows in §8. 

2. Equations cind numerical procedure 
2.1. Equations 

Our medium is an ideal compressible and stratified 
gas governed by the general equations of the magne- 
tohydrodynamics (MHD) including Ohmic diffusion: 



= -pV-v, 



Dp 
Dt 
Dv 

— = V X (v X B) r?AB 



(1) 



-Vp+ — (VxB)xB + pg, (2) 



(3) 
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-pV • V -t- (V X B)' , (4) 
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with e the internal energy per unit mass and rj the 
ohmic diffusivity, which is assumed constant. All 



other symbols have their customary meaning. Carte- 
sian coordinates {x,y,z) are adopted so that the z- 
direction is anti-parallel to g. In this paper we con- 
sider a two-dimensional problem: we assume that 
there are no variations of the physical variables along 
the y-axis [d/dy = 0), although By and Vy are gener- 
ally non zero. 

Equations (l)-(4) are cast in their conservative 
form and solved with a code written by Shibata 
1983. The latter is based on a modified Lax-Wendroff 
scheme (Rubin & Burstein 1967) and stabilized with 
artificial viscosity as described by Richtmyer & Mor- 
ton (1967). The inclusion of physical resistivity in our 
equations and the absence of any shock phenomena in 
the solutions permit us to minimize the use of the ar- 
tificial viscosity. The latter is restricted, in any case, 
to regions of very steep gradients only. 

This code has been repeatedly tested and used for 
two-dimensional simulations of the outbreak of the 
magnetic field at the surface of the Sun by Shibata 
and collaborators (see e.g. Shibata 1983; Shibata et 
al.l989; Kaisig et al.l990). In addition, we have suc- 
cessfully run several tests for our problem, checking 
for the convergence of the code as well as for the con- 
servation of mass, energy and magnetic flux in the 
box. 

In this paper, we only consider a case with left- 
right symmetry about a vertical plane containing the 
tube axis. For each half of the tube we use a numerical 
grid of 300 points in the horizontal direction and 700 
points in the vertical one (although in the figures pre- 
sented here, only a small fraction of the box is shown). 
In the following, the results are given in dimensionless 
form, using as units the background density, pressure 
scalehcight and Alfvcn speed calculated at the center 
of the tube at time t = 0. 

2.2. Initial conditions 

The initial condition consists of an unperturbed 
background atmosphere, with pressure and den- 
sity Pg, and, superimposed, a perturbation associated 
with a magnetic flux tube. To avoid any confusion 
due to pseudo-convective effects, the background at- 
mosphere is adiabatically stratified. It spans verti- 
cally 80% of the pressure scaleheight at the bottom 
of the box. The pressure contrast between the top 
and bottom is 2.6. The density contrast is 1.8. 

After inclusion of the magnetic tube, the result- 
ing system satisfies the following simple condition: 
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— VAp + J/cx B = 0, with Ap '= p — Pe the pressure 
excess as compared with the background stratifica- 
tion, J the electric current density and c the speed 
of hght. In the absence of gravity, this would be a 
perfect equilibrium condition. The density profile in 
the magnetic region is determined by assuming that 
the entropy in the tube is constant and equal to the 
unperturbed value in the atmosphere. The tube at 
time t = is thus buoyant [p < p^) and fulfills 



Ap def p- Pe 

P P 
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(5) 



to first order in with 7 the specific heat ratio. 
This case is intermediate between the two cixtreme 
possibilities of full thermal equilibrium, i.e., T = 
(which would be a factor 7 more buoyant) and the 
case of a tube with p — p^. The evolution pre- 
sented here is qualitatively very similar to the first 
case, whereas it deviates in important respects from 
the second. 

Along this paper we will deal with a number of 
different magnetic field profiles at time t = 0. The 
longitudinal field Bi = By will be taken to have a 
gaussian profile, 



Bi{t = 0) (X exp{ - / R^) . 



(6) 



The transverse field Bt{t = 0) is chosen to be purely 
azimuthal; calling r and <p the polar coordinates 
around the tube center, we have Bj(i = 0) = B^e^. 
For we will choose distributions such that the 
pitch angle 



*<^=^^atanf'^ 

Bi 



adopts an asymptotically flat proflle: 



B.ocB, 



ar" 



(7) 



(8) 



or, alternatively, an exponentially decaying one: 




The choice of n determines the rigidity of the tube 
center. In this paper we consider the case n = 3 
and a = 0.9. As we will see, both profiles (8) and 
(9) yield basically the same time evolution (§5.1). As 
explained in the discussion, we have also run some 
tests with fiat, top- hat magnetic profiles (§8.1). 



2.3. Boundary conditions 

The side and bottom boundaries are closed lids. A 
closed boundary at the top of the box must be avoided 
for two reasons: (1) the sound waves generated by the 
tube in the external medium must be able to leave the 
box and (2) the rise of the tube should not be unduly 
braked through the excess pressure at a closed bound- 
ary at the top. An ideal free boundary should be 
transparent for the outgoing disturbances and should 
not introduce disturbances through incoming waves. 
In the present calculations we achieve this by intro- 
ducing a fiducial layer (Nordlund & Stein 1990) well 
above the upper boundary of the box. At each time 
step, the pressure in the fiducial layer is calculated by 
assuming that each point in it is in static equilibrium 
with respect to the point on the boundary of the box 
lying directly below it; the densities are then obtained 
via the constant background entropy; the velocity is 
set to zero in the fiducial layer. The net mass flux 
across the boundary is not necessarily zero. However, 
the total mass of the box only varies by a maximum 
factor of 10~^ during each single run. 

3. The physical problem: parameters govern- 
ing the evolution 

The present physical problem is characterized by 
four basic dimensionless input parameters: (1) the 
thickness of the tube in units of the local pressure 

~ def 

scale height, R = R/Hp-, (2) the plasma beta at the 
center of the tube, /3o; (3) the ratio between the trans- 
verse and the longitudinal component of the field, as 
measured by, e.g., the pitch angle at a representative 
position of the tube, ^mt, to be specified in the fol- 
lowing; (4) the ohmic diffusivity rj in terms of, e.g., 

Vj^Hp, with Vj^ the Alfven speed: fj ri/{vj^Hp). 

These parameters determine the properties of the 
initial magnetic tube and are independent of the ve- 
locity of the flow that develops along time. However, 
in this problem there is a characteristic value for the 
flow velocity, namely the terminal speed of rise, v^^^^, 
given by the dynamical equilibrium between buoy- 
ancy and aerodynamic drag: 



^term 



7T R Ap 2 



(10) 



with Cd the customary drag coefficient, Co ~ 0(1); 
and Cg the sound speed. Hence, the customary 
velocity-related parameters, like, e.g., the Mach num- 
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bcr, the Reynolds number and the magnetic Weber 
number, can be immediately obtained as a function 
of the foregoing input parameters. For instance, using 
equation (5) we obtain that the Mach number must be 
of order: = 0{R/Po). As a result, we expect the 
rise to be very subsonic and, assuming ^ <C 1, sub- 
Alfvenic as well. This latter condition can be violated 
for tubes rising to higher levels of the convection zone. 

The rest of this section is devoted to a discussion 
of the values expected for the most important param- 
eters of this problem. 

3.1. The amount of twist and the deforma- 
tion of the tube 

The central parameter for the present paper is the 
twist of the tube as measured by the pitch angle of the 
field lines, 'I'. The values of this parameter of interest 
for our calculation are those for which the transverse 
component of the field, Bj, is able to suppress the 
conversion of the tube into a pair of vortex rolls. In 
the following we study different agents that tend to 
deform the tube. We come to the conclusion that the 
minimum value of Bf necessary to counteract them is 
the same for all. 

3.1.1. Differential buoyancy and pressure fluctua- 
tions along the boundary 

The tube of the initial condition explained in §2.2 
is increasingly buoyant toward its center. Thus, the 
central regions rise faster than the periphery; the up- 
per layers of the tube are thereby compressed whereas 
those located below the center are expanded. The 
magnetic field lines are deformed by this process; the 
transverse field, in particular, exerts an increasing re- 
sistance against further deformation. One can cal- 
culate the minimum transverse field that can effec- 
tively withstand (and reverse) the buoyant deforma- 
tion (Emonet & Moreno-Insertis 1996). At this point, 
it is sufficient to obtain an order of magnitude esti- 
mate for it that can serve as a guide for the rest of 
the paper. This can be easily done for tubes with 
smooth distributions like (6) and (8) or (9) and sat- 
isfying B? <C ■ The resulting criterion is best 
expressed in terms of the pitch angle calculated at 
the position of the maximum of the transverse field, 
*TOt; it reads: 



sin 



def 



> 



mt 



Ap 



2 



1/2 



For the tubes with initial condition §2.2, this ap- 
proximate criterion should hold within about a factor 
2. A more precise condition involving the detailed 
buoyancy distribution in the tube can be obtained 
by studying the magnetostatic equilibrium of a hor- 
izontal twisted tube with a non-homogeneous buoy- 
ancy distribution in the limit <C H^^ (Emonet 
& Moreno-Insertis 1996). It can be shown (see their 
Eqs. 44 and 49) that for the deformation of the tube 
to remain small, the radial profile of the pitch angle, 
^'(r), must satisfy: 



sin*(r) ^ 



< Ap{r) - (Ap) (r) ^ 



p{r) 



Pir) 



(12) 

where (/) (r) is the radial average of / between the 
radii and r. At each radius r, the deformation of the 
magnetic field is directly related to Ap{r) — {Ap) (r), 
i.e., to the differential buoyancy at that radius. 

The threshold due to the differential buoyancy 
(Eq. 11 or 12), however, cannot be the only criterion 
of interest for our problem. In fact, even a uniformly 
buoyant tube, if untwisted, is subjected to a defor- 
mation and conversion into vortex rolls. A deform- 
ing agent independent of the initial differential buoy- 
ancy is the pressure profile (i.e., the pressure fluctua- 
tions, Pg^, above the background stratification) built 
by the external flow around the boundary of the tube 
(Emonet & Moreno-Insertis 1996): a simple criterion 
for the resistance of the tube against those fluctua- 
tions is 



^1 



(13) 



Now, is itself of order the ram pressure of the 

external flow relative to the tube. Thus, it is no higher 
than about pi;|g^^/2. Substituting from Eq. (10), we 
obtain again a criterion as in (11). 

3.1.2. Vorticity generation 

Given the central role of the vorticity in the present 
problem, it is adequate to understand the minimum 
condition (11) in terms of the vorticity generation. 
The time evolution of the longitudinal component of 
the vorticity, cui, is governed by the following equa- 
tion: 



D I uji 



(11) 



Dt 
+ V 



X g 



V X ¥u 

P 



(14) 
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Fit represents the projection of the Lorentz force on 
the transverse plane {x,z). For definiteness, when 
speaking about the sign of the vorticity in the rest of 
this section, we refer to the right half of the tube only. 

At issue in this section is the generation of vortic- 
ity in the main body of the tube (the vorticity in the 
tube periphery is discussed in 5.2). The first term 
on the right hand-side is the counterpart in terms of 
vorticity of the gravitational torque: it produces posi- 
tive vorticity because V(Ap/p) points outwards. The 
second term on the right-hand side of (14) represents 
the effect of and can be rewritten as: 



VxFi, _ (B,.V)J; 



cp 



(15) 



with J; the longitudinal component of the electric cur- 
rent density. At the beginning of the run, the longitu- 
dinal current is axisymmetric so that V x Fit is zero. 
As soon as the tube center begins to rise relative to 
the periphery, the longitudinal current is enhanced in 
the upper half and diminishes in the lower half of the 
tube. According to Eq. (15), this produces negative 
vorticity in the interior of the tube, i.e., it tends to 
counteract the effect of the gravitational torque. By 
setting these two terms equal, we obtain a criterion 
for the minimum transverse magnetic field that can 
effectively oppose the initial deformation of the tube: 
the result is, again, (11). For the problem we are con- 
sidering here, the last term of the vorticity equation 
(14) is 0{Ap/ p + L/jHp) smaller than the other two 
terms and is therefore not of primary interest here {L 
represents a local characteristic length for the trans- 
verse magnetic field , e.g., the radius of the tube or 
the thickness of the tube boundary). 

Criteria equivalent to (11) can also be obtained 
through other physical considerations. For instance, 
Tsinganos (1980) obtained a similar threshold for \1/ 
by calculating the stability of the tube against split- 
ting due to the development of Rayleigh- Taylor and 
Kelvin-Helmholtz instabilities at its apex. 

3.2. Other parameters 

3.2.1. Plasma beta and tube radius 

The choice of values for the parameters /3o and R 
must be guided both by (astro)physical insight and 
feasibility of the numerical calculation. /?o is expected 
to be very high in the solar interior, like, e.g., O(IO^) 
for the magnetic tubes at the bottom of the convec- 
tion zone. The Courant condition for the numerical 



code, though, sets a stringent upper limit to the value 
of Po that we can use. Using expression (10) we can 
calculate the number of timesteps necessary for our 
tube to cross the whole integration box, nt, as 



1/2 



(16) 



with rig the total number of cells in the vertical direc- 
tion. A trade-off between field intensity and magnetic 
flux of the tube is then necessary. As a compromise 
we have chosen /?□ = 10^ and R = 3.7 10"^. This 
should be sufficient to understand qualitatively the 
flows in and around the tubes with the fleld strength 
and flux expected in the progenitors of active regions. 

3.2.2. The Magnetic Reynolds number, Re^ 

The magnetic Reynolds number can be defined in 
this problem in terms of the terminal speed of rise 
and the tube radius. Its order of magnitude is thus 
given by Re^ = 0{R^/^ /fj). The choice of a value for 
this parameter must be guided primarily by numeri- 
cal criteria: large magnetic field gradients develop in 
the tube along the evolution, specially at its upper 
rim. We have chosen the ohmic diffusivity so that 
the resulting value of i?e„ is a few times 10^. This 
yields good numerical performance but causes some 
unwanted diffusion of the field outward from the tube. 
A measure for the latter can be obtained from the y- 
component of Walen's equation, obtained by combin- 
ing the continuity and induction equations, (1) and 
(3): 

In the central region of the tube, the variation in Bi/p 
is basically due to diffusion. The calculations pre- 
sented in the following are not diffusion-dominated 
(as expected from the high Re^), in the sense that 
the decrease of this quantity is small in the tube cen- 
ter during the period of time shown in the figures. 
This can also be checked through an order of mag- 
nitude estimate: the diffusion term in (17) for the 
initial gaussian distribution can be compared to the 
rate of change of the field (left-hand side of the same 
equation) imposed by the rise of the tube. The ratio 
between them is {Re^Ry^, which is 0(10"^) in the 
present case. 

3.2.3. The Weber number, We 

In the hydro dynamical literature on air bubbles. 
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Fig. 1. — Rise of a twisted tube through a stratified environment. The initial piteh angle "^nu is 7°. The grey scale 
and the contours correspond to the longitudinal magnetic field intensity. At each time step white corresponds to 
the maximum of Bi at this instant and black to 1 percent of this maximum. The arrows correspond to the velocity 
field. In the figure, only a about 1/5 of the total integration box is reproduced; more precisely, 100x500 out of a 
total of 300 X 700 points for each half of the figure are included. The times represented are: t = 0, 1.3, 2.5, 3.8, 5.3, 
6.6. The complicated shape of the wake in the two last panels is due to the episodic release of small secondary rolls 



the Wcbor number is used to measure the relative 
importance of the inertial forces of the flow to the 
surface tension at the boundary of the bubble (e.g. 
Ryskin & Leal 1984b). The role of surface tension is 
played in our case by the jump in magnetic tension 
of the transverse component at the boundary of the 
tube. Hence we define the magnetic Weber number 
as 

def 



We 



BtV(47r) 



(18) 



Inserting for v the terminal speed we easily obtain 
We = .R/3o |Ap/p|o/(2sin^ *mt)- The condition on 
the minimum pitch angle to avoid the splitting of the 
tube, (11), can thus be reformulated as the condi- 
tion that the Weber number be at most of order one, 
We^l. 

4. The rise of a tw^isted tube in a moder- 
ately stratified environment: a represen- 
tative case, = 7° 

The critical value of the pitch angle parameter '^mt 
obtained by substituting in (11) the chosen values of 
/3o, R and Ap/p is 6°. In the present section we 
describe some major features of the rise of a tube with 
a pitch angle "ifmt close to that value, ^mt = 7°. The 
time evolution is illustrated in Fig. 1, which shows a 
fraction of the integration box used (30% and 70% 
of the box in the horizontal and vertical directions, 
respectively). In §7 we compare the results for flux 
tubes with different values of the initial '^mt- 

4.1. The initial acceleration phase 

In the absence of external flows to keep it in place, 

the tube lacks equilibrium globally and starts rising, 
at the same time sending out sound waves all across 
the box. The initial global acceleration is basically 
free-fall, i.e., (g/i) (Ap/p), with an overbar indicat- 
ing average values in the tube and i standing for the 
enhanced inertia factor, which turns out to be 2 (see 
§4.3.1). 

Superimposed on the global rise, different kinds of 
motions within and around the tube take place which 
tend to deform its initial axisymmetric shape (Fig. 1, 
upper row) . One of them is the faster rise of the tube 
center compared with the periphery due to the dif- 
ferential buoyancy (§3.1.1). As a consequence of this, 
the magnetic tension associated with the transverse 
fleld builds up at the tube front; the relative motion 
of the tube center is thus stopped and reversed. A 



Fig. 2. — Distance to the stagnation point at the tube 
front of 50 mass elements initially located at equidis- 
tant points along the vertical symmetry axis as a func- 
tion of time(solid lines). The compression of the tube 
front and the oscillations of the central regions along 
the vertical axis can be clearly seen. The curves plot- 
ted vertically (stars) are profiles of the transverse field 
intensity, Bt, along the vertical symmetry axis at reg- 
ular intervals of time. 



vertical oscillation of the tube center ensues within 
the tube's cross section. The compression/expansion 
in the upper/lower half of the tube together with the 
internal oscillations are clearly visible in Fig. 2, where 
we have plotted the position as a function of time of 
several Lagrange markers located along the vertical 
symmetry axis of the tube. 

Simultaneously with the processes just explained, 
the external matter slides around the tube and drags 
toward the rear the outermost tube layers where the 
magnetic field is too weak to oppose any important 
Lorentz force. The field is thus stretched all along 
the boundary (§5). Vorticity is being generated in 
the matter being dragged: as a result, two magne- 
tized vortex rolls are created that trail the tube mo- 
tion (Fig. 1, third panel). The main body of the tube 
thereby loses about 30% of its original magnetic flux. 
This figure sensitively depends on the initial pitch an- 
gle (§7). Given the sign of the vorticity, the matter 
between the rolls is moving upwards with respect to 
the back of the tube. A pressure excess appears at 
the upper end of the inter-vortex space, directly be- 
low the lowermost tube layers. 
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Fig. 3. — Close-up view of the tube interior show- 
ing the velocity vectors of the internal torsional os- 
cillation. The arrows in the figure correspond to the 
relative velocity of the individual mass elements with 
respect to the tube apex. 



4.2. Torsional oscillations of the tube interior 

The vertical oscillations of the tube center are in 
fact part of a torsional oscillation in which most of 
the tube interior is taking part (Fig. 3). The torsional 
oscillation has left-right symmetry: each tube half is 
rotating back and forth around a horizontal axis offset 
by a fraction of R from the midplane. The energy of 
this oscillation is being radiated away from the tube 
via pressure forces; it is also damped through the dif- 
fusion (physical and numerical) present in the code. 
Thus it slowly decreases in amplitude. The frequency 
of these oscillations is of order ivtors = 0[vAt/ (27rii!)]. 
Hence, one expects a few torsional oscillations to be 
completed while the tube is rising across a scaleheight. 

4.3. The asymptotic regime of rise 

The strong initial acceleration phase is followed by 
a quasi-stationary asymptotic regime (Fig. 1, lower 
row). In it, the rate of change of the tube interior 
and external flows becomes small compared with the 
initial phases. The transition occurs at about t = 3 
(between the third and fourth panels of Fig. 1). The 
asymptotic regime is characterized by the adoption of 
a terminal speed of rise, a sharp separation of tube in- 
terior and surroundings and a well-developed trailing 
wake. Simultaneously, the tube stops losing magnetic 
flux to the wake and there is no further important 
deformation of the magnetic fleld lines in the head of 
the tube. The following subsections deal with some 
of the features characterizing this asymptotic regime. 



4-3.1. The terminal velocity 

Once the wake is developed, the resistance of the 
surroundings to the advance of the tube can be cal- 
culated using the customary aerodynamic drag force 
with almost constant drag coefficient Cd- In the 
asymptotic regime, the total buoyancy of the tube 
is a slowly varying quantity. Thus we can expect the 
tube to reach a terminal speed of rise given by: 



7) — 



2gSAp 

Cjjdpe 



1/2 



(19) 



S and d being the cross sectional area and horizon- 
tal diameter, respectively. The order of magnitude of 
this speed was estimated in Eq. (10). The best flt to 
the numerical results yields a value of Cjj of approxi- 
mately 1.6. The approach to this asymptotic regime is 
not immediate to formulate analytically, since the co- 
efficient Cd is not constant while the tube is changing 
its shape and developing a wake. The initial acceler- 
ation of the tube as a whole is 



(20) 



the factor i allows for the added inertia due to the co- 
acceleration of the external medium (i = 2 for a rigid 
straight cylinder). To approximate the velocity of the 
tube in the initial and intermediate time-dependent 
phases, we use the expression for a tube which starts 
from rest with acceleration a, is being acted upon by 
a constant driving force and has fixed Cd, viz.: 
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Fig. 4. — Velocity of the tube center (stars) and of 
the apex (circles) as a function of time. The dot- 
ted line corresponds to Eq. (21). The dashed line 
represents the initial acceleration of the tube center, 
{\Ap/p\)r=o g/i, with i = 2. 
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a constant driving force and has fixed Cd, viz.: 



tanh 



\ ^term / 



(21) 



with t^term ^tiH given by Eq. (19). 

The numerical results show a surprising closeness 
to Eq. (21). In Fig. 4, we show with a dotted line the 
velocity given by Eq. (21), while the actual speed of 
the tube apex, Wapg^., is indicated as circles. As can 
be seen, the tube speed oscillates around the mean 
speed given by Eq. (21). This is due to the strong os- 
cillations of the tube interior: in fact, the tube center 
oscillates with a much larger amplitude, as is shown 
by the curve with asterisks. The tube center pushes 
the apex, so that circles and asterisks arc slightly out 
of phase. The amplitude of the oscillation diminishes 
as its energy is being radiated to the external medium 
and both curves converge toward the dotted line. The 
initial acceleration of the tube apex is very close to 
(20) with i = 2, i.e., as corresponds to the global buo- 



yancy of the tube. The tube center, in turn, has a 
higher acceleration, corresponding to the local value 
of the buoyancy, namely {\Ap/ p\)r=o g/i (this value 
is shown in the figure as a dashed line). Here again, 
the enhanced inertia factor i must be set equal 2. 

As a test to the validity of the terminal velocity for- 
mula, we have calculated the evolution of tubes start- 
ing with the same initial condition as the one used to 
obtain Fig. 4 (in particular, the same tube radius, R) 
but with Pq chosen in a range between 1000 and 100. 
The velocity curves of all these tubes are similar to 
those shown in Fig. 4, with the time scale contracted 
by a factor /Jg^^ and the velocity scale enlarged by 
the same factor. More precisely, if all these tubes 
reach the aerodynamic drag regime with terminal ve- 
locity given by Eq. (19), then the curves depicting the 
velocity of the apex normalized by the plasma beta, 

1 /2 

v^p^^Pq , versus the position of the apex must all 
basically coincide (the drag coefficient Cjj does not 
change much with varying Reynolds number in this 
regime). The results of the tests show indeed a large 
degree of superposition of those curves: the maximum 
relative deviation between them is a factor 1.1, which 
can be attributed to the change in Cd consequent 
with the variation in the Reynolds number. 



Fig. 5. — Equipartition line a.t t = 3.8 together with 
the flow field relative to the apex of the tube. The 
simple structure of the wake at this time is clearly 
visible. The envelope of tube and wake has a well 
defined elliptical shape. This figure bears a strong 
resemblance to the experimental results of air bubbles 
rising in liquids (e.g., Collins 1965, also reproduced by 
Batchelor 1967, plate 15). 



4-3.2. The equipartition line at the tube periphery 

As a result of the resistance to deformation pro- 
vided by the transverse field, there is a neat sepa- 
ration between the tube interior and the outside in 
the asymptotic regime. A good indicator of the loca- 
tion of the tube boundary can be obtained by com- 
paring the kinetic energy density of the relative flow 
field, v^^; = V 



• fifein = with 

the energy density of the transverse magnetic field. 



apex ' 



dcf 



= Bi /(Stt). We define the egMzparfziicm, Z/'ne as 
the locus of those points in the tube periphery where 
the two energy densities are equal. Fig. 5 shows, at 
time t = 3.8, the relative flow field v^^^ and, super- 
imposed on it, the equipartition line. The flow field 
changes markedly when going across the equiparti- 
tion line: the velocities inside are much smaller than 
outside it. In other words, once the terminal velocity 
has been reached the tube rises basically as a unity, 
with only a weak internal flow pattern corresponding 
mainly to the torsional oscillations described in the 
previous section. A comparison of Fig. 5 with the 
results of laboratory experiments (e.g., Collins 1965 
or plate 15 in Batchelor 1967), shows that there are 
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Fig. 6. — Profile of the total pressure excess along 
the boTinciary of the tube, which is defined by the 
equipartition line (see §4.3.2), at time t = 2.5. The 
boundary of the tube (the equipartition line) is shown 
clS cL thick line in the left panel. The feature in the 
middle of the profile (between 110° and 130°) results 
from the shape of the equipartition line at the lat- 
eral edge of the tube. The arrows on the left panel 
correspond to v^^^ and the radial lines are drawn at 
constant azimuthal angles. 

striking similarities between buoyant magnetic tubes 
and air bubbles rising in a liquid, both in the shape 
of the rising object itself as in the wake, in spite of 
the very different parameter values involved (like, e.g., 
the density deficit). 

All along the initial phase, a pressure profile is set 
up around the tube periphery which closely resembles 
the pressure distribution around rigid cylinders in rel- 
ative motion with the surroundings. To show this, we 
can plot (Fig. 6) the pressure in the points along the 
equipartition line. In this figure, the pressure excess 
associated with the return flow between the vortex 
rolls of the wake is clearly visible. 

5. The sharp interface at the tube boundary: 
magnetic boundary layer 

All around the equipartition line there is a region 

with non-zero vorticity. This contains (1) a bound- 
ary layer surrounding the tube and (2) the trailing 
wake. Further out, there is the external medium with 
a largely vorticity-free flow. In this section we study 
the structure of the boundary layer; the wake is con- 
sidered in §6. 

The two main physical features characterizing the 
magnetic boundary layer around the tube are a strong 
shear of the tangential flow and high magnetic field 
gradients. They lead to enhanced generation of vor- 
ticity, ohmic diffusion and generation of entropy. In 



the following subsections we study these items in turn: 

5.1. Shear flow and field gradients 

A pronounced shear is taking place in a thin band 

around the equipartition line (Fig. 5). As a result, 
the matter elements and the transverse field are be- 
ing stretched at the tube periphery all around the 
tube. To visualize this phenomenon, we use Lagrange 
markers which follow the motion of the tube's mass 
elements (Fig. 7). We choose six groups of markers at 
time t = Q (left panel), at a distance of the tube cen- 
ter such that Bi is 5% (white asterisks), 27% (black 
triangles) and 77% (black diamonds), respectively, of 
-B/Q (the maximum value of Bi in the tube) both in 
the upper and lower half of the tube. The markers 
are set in the neighborhood of the tube axis and are 
left to evolve with the tube. As time advances, the 
markers of the outermost group in the upper tube half 
are stretched by a very large factor along the tube pe- 
riphery (center and right panel of the figure); in fact, 
many of them are brought all the way down to the 
wake. The markers of the group immediately below 
them (triangles) are brought close to the uppermost 
group and stretched, albeit by a smaller factor than 
the other group. Large compression and stretching 
happen basically only in the neighborhood of the in- 
terface between tube and surroundings: closer to the 
tube center (diamonds) the mass elements are just 
periodically moved following the internal torsional os- 
cillation. At the rear of the tube, the Lagrange mark- 
ers with Bi = 0.275(0 (triangles) are also somewhat 
stretched by the flow of the trailing wake (see §6). 
The outermost parcels at the back of the tube (white 
asterisks) are trapped in the tail and no longer rise 
with the tube. 

Corresponding to this shearing, the pitch angle in- 
creases to large values at the tube periphery: in Fig. 8 
we plot the pitch angle distribution along the verti- 
cal axis of symmetry for the three instants shown in 
Fig. 7. The pitch angle becomes large in the inter- 
face at the tube front (middle panel). In the quasi- 
stationary regime, its distribution shows a horizontal 
asymptote outside the tube similar to the one of the 
initial condition. This feature is not a mere conse- 
quence of our initial condition (8). In fact, an initial 
pitch angle distribution with exponential decline for 
large radius at time t = (see Eq. 9) yields in the 
quasi-stationary regime a horizontal asymptote sim- 
ilar to the one of Fig. 8. It is then perhaps more 
natural to choose the initial condition of Eq. (8) from 
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Fig. 7. — Motion of individual mass elements in the tube, showing the stretching of the matter and field along the 
tube boundary. The panels correspond to t = 0, 1.9 and 3.8. At time zero the Lagrange markers are located at 
R/Hp = 0.044 and 0.068 where Bi is equal to 27% and 5% of its value at the center. Radial compression, azimuthal 
stretching and shearing of the upper layers of the tube is apparent. The elements indicated by black triangles at 
the back are also somewhat stretched as a result of flow in the wake. 



the outset. 

The large pitch angles apparent on the right of 
Fig. 8 are all in a region of the tube where the mag- 
netic field intensity is very low; the central regions, 
on the other hand, have low pitch angles. Thus, one 
should not expect a global kink-unstable behavior of 
the tube. This can be checked by calculating the 
net tension in the tube over background (see §9.2 in 
Parker 1979). With the convention that tension is 
positive and pressure negative, we obtain a positive 
value indicating that the tube is under longitudinal 
tension rather than compression. Thus, there should 
be no tendency to buckling of the tube. 

5.2. Vorticity generation in the boundary layer 

The large jump in tangential speed visible across 
the equipartition line in Fig. 5 marks the presence of 
a vortex sheet surrounding the tube. In fact, vort- 
icity (more precisely, its longitudinal component uj]) 
is generated at the boundary layer during the whole 



duration of the run. From there, it is advected toward 
the wake. 

These processes are governed by Eq. (14). As in 
§3.1, the right-hand side is dominated by the first 
two terms [the third is again 0(l//3o) smaller]. How- 
ever, along the front of the tube those two terms, i.e., 
the gravitational and magnetic vorticity sources, rein- 
force each other rather than mutually cancel, as in the 
tube interior. This is because the radial derivative of 
the transverse field component reverses its sign in the 
neighborhood of the equipartition line. For instance, 
in the simple case in which the field lines are still not 
far from circular (e.g., along the front) the magnetic 
contribution is approximately, 



p cp d(j) 



r or 



(22) 
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Fig. 8. — Pitch angle distribution (stars) along the 
vertical central axis at three different times corre- 
sponding to the 3 panels of Fig. 7. The profiles of 
Bi (dashed) and Bt (dotted-dashed) are overplotted. 



so that J; changes sign at about the maximum of Bt. 
As a result, (and taking into account the sign of the 
azimuthal derivative d/dcj)), in the right half of the 
tube positive vorticity is generated all along the front 
of the tube, and negative vorticity is generated at the 
rear, i.e., at the interface between the tube and the 
wake. At the rear \/{Ap/ p) is approximately parallel 
to g and therefore the magnetic stress is anyway the 
main source of vorticity there. 

In Fig. 9 we show the distribution of vorticity in 
the box at four different times during the evolution. 
White and black indicate positive and negative vor- 
ticity, respectively (i.e., clockwise and anticlockwise 
rotation), whereas the grey background is the zero 
vorticity level. Soon after the beginning (leftmost 
panel; t = 0.3) the gravitational torque; is creating 
positive and negative vorticity in the right and left 
tube halves, respectively, across the whole tube inte- 
rior. In the second panel {t = 1.6), this has already 
been countered by the magnetic tension of the trans- 
verse field: the reversal of the sign of the vorticity 
in the interior reveals the first backward torsional os- 
cillation. Additionally, a vortex sheet is already ap- 
parent all around the tube front. At the rear, there 
is only an incipient wake and no clear vortex sheet 
yet. The third and fourth panels correspond to the 
asymptotic phase (t ~ 3.8, i.e., close to the beginning 
of the stationary regime, and at an advanced stage, 
t = 5, respectively). Vorticity generation now occurs 
basically close to the equipartition line only. Vortex 
sheets along both the front and the rear of the head 
of the tube are visible. As expected, the sign of the 
vorticity in them is opposite. A detailed distribution 
of oji in the boundary layer around the tube (from top 
to bottom), can be seen in Fig. 10. The vorticity in- 
deed changes sign at the edge of the tube: the abrupt 
shape of the zero crossing is due to the acute shape 
of the latter. The small maximum at the right end of 
the distribution is characteristic of the tail (§6.2). 

5.3. Ohmic diffusion and field advection 

The structure of the magnetic boundary layer is 
dominated by the ohmic diffusivity, which in part 
plays in the present case the role of the viscosity in 
the standard hydrodynamic boundary layers. An im- 
portant difference between viscous and purely mag- 
netic boundary layers, though, concerns the no-slip 
condition for the flow: in the purely magnetic case, 
there is no counterpart to the viscous second-order 
derivative term that could force the tangential veloc- 
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Fig. 9. — Distribution of the y-componcnt of the vorticity vector in and around the tube at t = 0.3, 1.6, 3.8 and 
5.0. White corresponds to positive vorticity (clockwise rotation) and black to the opposite. The two leftmost panels 
clearly show the vorticity associated with the torsional oscillations. The boundary layer and vortex sheet at the 
tube front arc visible in the three rightmost panels. In the two last panels the vortex sheet at the interface between 
the wake and the tube is also visible. 



ity to reach a zero value at the boundary. In our case, 
however, the flows on the inner side of the equiparti- 
tion line are reduced or suppressed by the magnetic 
forces. In other words, in the present boundary layer 
the Lorontz force imposes a no-slip condition which 
does not allow a penetration of the external incoming 
flow into the tube nor large fluctuations of the trans- 
verse velocity inside. This is clearly visible in Fig. 11 
where we have plotted v^, {v^gi)z, -Bt and Bi along a 
vertical axis slightly offset from the symmetry axis. 

The slope of the profile of Bt and Bi across the 
boundary layer results from the equilibrium between 
the outward ohmic diffusion and the field advection 
through the external flow. For Bt, for instance, this 
is controlled by the transverse component of the in- 
duction equation (Eq. 3). Written in terms of the 
magnetic potential A, the latter reads: 

{v^^i • V) A = AA , with Bt = V X (Aiy) . 

(23) 

As the transverse field diffuses outward, it swallows 




Fig. 10.- Profile of the longitudinal component of 
the vorticity (right panel) along the cquipartition line 
drawn in the left panel for t — 3.8. The zero cross- 
ing at about 120° corresponds to the transition from 
tube front to rear at the edge of the tube. The radial 
lines in the left panel are drawn at constant azimuthal 
angles. 

some of the incoming external matter (which, in this 
way, slowly enters the magnetic field system). Simul- 
taneously, this same external matter advects the mag- 
netic field back toward the tube interior. The equilib- 
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Fig. 11. — Plots of Bi and (dash-dotted) together 
with Uj. (solid line) and {Vj.^i)z (dashed) along a ver- 
tical axis a,t X = 0.01. The vertical dotted lines cor- 
respond to the positions where 
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rium between those two terms is reached for a given 
thickness of the boundary layer, ij,;. To obtain an 
estimate for Lj,;, we note that the external flow has a 
lengthscale comparable with the tube radius (a com- 
mon feature of cylinders moving in fluid media), so 
that, 

R 



(24) 



This estimate is in agreement within a factor about 
2 with the results of the numerical calculation (see 
Fig. 11). An analogous advection-diffusion equilib- 
rium can be seen to hold for Bi in the boundary layer. 



5.4. Comparison with air bubbles and rigid 
tubes. 

In spite of the striking similarities between the 
present problem and the rise of an air bubble, the 
corresponding boiuidary layers at the periphery of the 
rising object differ in important ways. The boundary 
of an air bubble is basically a free surface for the ex- 
ternal medium, since the density and viscosity in the 
interior are negligible compared with those outside. 
The external flow must then satisfy a zero tangen- 
tial stress boundary condition rather than the no-slip 
condition of the present magnetic tube (e.g., Batch- 
elor 1967, §5.14; Ryskin & Leal 1984a, 1984b). As a 
consequence, the tangential velocity at the gas/liquid 



interface in the air bubble docs not vanish, but must 
instead satisfy the following equation: 



(25) 



with K the curvature of the boundary layer and 
(^rei)ts tangential component of the relative flow 
speed. The resulting jump in tangential velocity 
across the boundary layer of a bubble is therefore 
only of order Re~^^^, while it is 0(1) for a no-slip 
boundary layer, as in the present problem. 

The magnetized boundary layer of the rising tube, 

in fact, resembles more closely the boundary layers 
around solid bodies with a no-slip condition than 
those around air bubbles. The distribution of vor- 
ticity along the equipartition line shown in Fig. 10 is 
quite similar to that around the boundary of a rigid 
cylinder at Re ^ 300 (e.g., Ta Phuoc Loc 1980). In 
our case, though, the reversal of the sign occurs more 
abruptly because of the particular geometry of the 
equipartition line. 

5.5. Entropy generation 

The enhanced ohmic dissipation in the boundary 
layer is creating entropy at a rate: 

Dt\Cp J 7 47rp ' ^ ' 



To order of magnitude, this entropy increase can be 
written: 

^ f^A ^ 1 '"rise (07\ 

Dt[c,)^/3,, R ' ^^'> 

with the local plasma beta in the boundary layer. 
The entropy increase for a matter element moving 
all along the tube boundary is then of order 1//?^;. 
As a result, the density difference Ap/p of those ele- 
ments can be substantially modified. This, however, 
is unlikely to cause important modifications in the 
dynamical behavior of the tube boundary. The main 
driving force there, the pressure gradient, is of or- 
der pw^jgg/i?; therefore, it should be more important 
than the local gravitational term gAp by the ratio 
Bq/B^I » 1, where is the magnetic field strength 
in the boundary layer. 

In our numerical simulation we also have some in- 
crease of entropy due to the artificial viscosity. Its 
infiuence, however is only secondary because of the 
low value of Re^ chosen. 
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6. The wake 

Wakes have been the object of active research in 
fluid dynamics for at least the past fonr decades. This 
includes the wake behind solid cylinders (e.g., Collins 
& Dennis 1973a; CoUins & Dennis 1973b: Ta Phuoc 
Loe 1980, Bouard & Coutanceau 1980; Ta Phuoc Loc 
& Bouard 1985), drops (Dandy & Leal 1989; Stone 
1994) and air bubbles (Davies & Taylor 1950; CoUins 
1965; Parlange 1969; Wegener & Parlange 1973; Hnat 
& Buckmaster 1976; Ryskin & Leal 1984b; Christov 
& Volkov 1985). In this chapter we briefly describe 
the formation and structure of the wake behind the 
rising magnetic tube and pay special attention to the 
similarities and differences to other wakes known in 
the literature. 

6.1. The time evolution of the wake 

Solid cylinders, air bubbles and drops have a clearly 
defined boundary from the beginning. Vorticity is cre- 
ated along the boundary and advected downflow. The 
wake is formed out of external fluid only, via, e.g. in 
solid bodies, the detachment of the boundary layer. 
In our case, in contrast, there is no such clear bound- 
ary at the beginning. The layers of weak magnetic 
field in the outskirts of the tube (a) are unable to re- 
sist the incipient external flow so that they are bodily 
convected to the rear and (b) increase their vortic- 
ity through the mechanisms explained in §3.1.2 and 
§5.2. It is the accumulation of this rotating material 
at the back of the tube that constitutes the wake. The 
material making up the wake may in fact be a large 
fraction of the initial tube in the case of small initial 
pitch angles (§7). 

The process of formation of the wake is basically 
complete by t = 3.0 (sec Figs. 1 and 5) and this sig- 
nals the transition to the asymptotic regime of rise 
(§4.3). The latter, however, is not exactly station- 
ary, since the wake does not remain unmodified in 
the rest of the evolution. In fact, it slowly grows in 
size and episodically changes its shape. The growth 
can be seen in Fig. 12, where the vertical size of the 
rolls is depicted as a function of time. The elongation 
is roughly linear, with a small-amplitude oscillation 
superimposed. This kind of behavior is not unknown 
in fluid dynamics. The elongation of the wake is in- 
deed similar to what is observed behind a rigid cylin- 
der moving at constant speed with = 300 after 
an impulsive start (see Ta Phuoc Loc 1980): in the 
laboratory experiment the linear growth of the wake 
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Fig. 12. — Vertical size of the wake as a function of 
time, calculated as the distance between the stagna- 
tion point at the rear of the tube and the stagnation 
point at the bottom of the wake. The growth is quasi- 
linear, with oscillations superimposed. 



with time continues well after the high-i?e flows have 
been set up around the cylinder and Cd has reached 
a constant value close to 1. For higher values of Re 
{— 3000) the growth of the wake with time is expo- 
nential rather than linear (Ta Phuoc Loc & Bouard 
1985). In our case, a small part of the growth of the 
wake is due to the stratiflcation; however, its main 
cause is the transport of rotating material from the 
boundary layer along the boundary of the tube to the 
wake. This is clearly visible in the third and fourth 
panels of Fig. 9: a white streamer of positive vortic- 
ity coming from the tube front is being wrapped up 
first around and then into the roll on the right-hand 
side of the wake. Simultaneously, a black tongue with 
origin in the back of the tube is also making inroads 
into the wake, next and following a parallel path to 
the white streamer. 

In the case of rigid tubes and bubbles moving 
through a fluid (and assuming perfect symmetry about 
the object's midplane), the growth of the wake nor- 
mally lasts until a dynamical equilibrium is adopted. 
In it, the vorticity generated in the boundary layer 
is advected to the wake, where it is redistributed 
through viscous diffusion and further advection. Si- 
multaneously, the energy gained through the driv- 
ing force (the potential energy for the bubbles) is 
transformed into heat in the wake, again through vis- 
cous dissipation (Parlange 1969; Wegener & Parlange 
1973). Further changes may occur if there is no mech- 
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Fig. 13. — Formation and release of a secondary roll at the edge of the tube. The panels contain a time series 
showing the equipartition line and the relative flow field. The secondary roll is being released to the downflow; later 
on, it gets dragged all around the wake and noticeably perturbs the shape and velocity distribution of the latter. 
Note also the jumps in velocity inside the main roll of the wake, revealing the presence of vortex sheets which have 
been advected from the tube front. 



anism to dispose of the heat thus generated. 

In the results presented here, no dynamical equilib- 
rium is reached since there is no equivalent mechanism 
to diffuse the vorticity effectively. This can be seen 
quantitatively as follows: the characteristic time for 
the adoption of a stationary state in the wake through 
viscous diffusion and dissipation is Re {L^g^^^/ Rf' 
times longer than the characteristic time for the for- 
mation of the wake itself, L^^^^ being the transverse 
dimension of the wake and Re the viscous Reynolds 
number. In our case, L^^f.^ « R and we expect Re 
(e.g., calculated on the basis of the numerical viscos- 
ity) to be at least as high as -Re„. Therefore, the es- 
tablishment of the dynamical equilibrium lasts much 
longer than the evolution presented in this paper. The 
consequence of this is the very inhomogencous vor- 
ticity distribution apparent in Fig. 9. In a high-i?e 
case, in fact, it may never come to a full dynami- 
cal equilibrium because an instability interrupts the 
approach to a steady state: the two standing rolls 
develop asymmetric oscillations until one of them is 
released downstream (see, e.g., Batchclor 1967, plate 
10). In contrast, in low-i?,e bubble experiments, the 
trailing vortex rolls are similar to a Hill's vortex (e.g., 
Colhns 1965; Parlange 1969; Wegener & Parlange 



1973; Ryskin & Leal 1984b; Christov & Volkov 1985). 

6.2. Further structural features common with 
laboratory experiments. 

6.2.1. The skirt 

In experiments with buoyant air bubbles, it is of- 
ten observed that a small quantity of air is torn 
off (or sucked from) the bubble corner and becomes 
stretched all around the wake. The result, commonly 
called the skirt, is a sheet of air aligned with the in- 
terface between the wake and the external medium, 
which remains attached to the bubble (e.g., Hnat & 
Buckmaster 1976). A similar feature can be identi- 
fied in our numerical simulations coinciding with the 
streamers or tongues of vorticity of both signs men- 
tioned in the last section. In contrast to the air bub- 
bles, the skirt here is made of magnetized material. 
During the whole rise the magnetic field in the skirt 
is continuously compressed and stressed both by the 
wake and external flows. 
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Fig. 14. — Tube and wake structure for four cases with different initial pitch angle but otherwise equal initial 
conditions. From left to right, the values of '^mt were 13.9°, 7° (i.e., the case studied in the rest of the paper), 
2.5° and 0°. The lower panels show the general flow field and the contours of magnetic intensity for S;. The upper 
panels represent the velocity relative to the apex of the tube and the corresponding equipartition line (as defined in 
§4.3.2). The transition from a convex tube back to a flat and then to a concave one for decreasing initial "i^mt (i-S-, 
as the Weber number increases) is apparent. In the untwisted case, no head is left but, rather, a bridge joining the 



6.2.2. The tail 

All along the central symmetry plane, the rising 
tube is leaving behind a tail of weakly magnetized 
matter (Fig. 1). Its upper part is compressed by the 
two rolls of the wake. In fact, the tail is not only 
magnetized, but also has non-zero vorticity. This can 
be seen both in the distribution of Fig. 10 (maximum 
of positive vorticity at the right end of the figure) 
and, as a light shadow, in the vorticity map of Fig. 9. 
Similar structures are seen behind moving pair of vor- 
tices as well as in bubble experiments (e.g. Wegener 
& Parlange 1973). The total amount of magnetic flux 
which is lost by the tube to the tail is small (only 6 
percent of the total flux). 

6.2.3. The secondary rolls at the edge of the tube 

Close to the lateral edges of the tube (see Fig. 13 
and the last panel in Fig. 9) small secondary rolls are 
built from time to time; they grow for a while and 
then they are released to the downflow. These fea- 
tures are reminiscent of the secondary rolls observed 
near the detachment point of boundary layer around 
solid cylinders (Ta Phuoc Loc 1980; Ta Phuoc Loc & 
Bouard 1985). These secondary rolls strongly perturb 
the skirt and the wake when they are released: this is 
the main cause for the irregular shape of the wake in 
the last panel of Figs. 1 and 9. 

7. The rise of tubes with different initial de- 
grees of twist 

The initial value of "^mt for the tube whose rise we 
have described in the foregoing chapters was close to 
the value given by the criterion (11), namely '^mt = 
6°. This criterion provides a lower bound for the 
twist necessary to withstand the destruction of the 
tube through vorticity generation in its interior, ex- 
ternal pressure forces, etc; we do not know yet how 
sharp that bound is. In the present section we inves- 
tigate the transition from the low-twist to the high- 
twist regimes. We conclude that the change from one 
to the other is gradual and takes place over a range 
of, say, 5 degrees in '^mt, around the value given by 
the right hand side of Eq. 11. In Fig. 14 we show the 
relative flow, v^^; (upper panels) and the contours 
of longitudinal magnetic intensity, Bi (lower panels), 
for four cases with = 13.9° ,7° ,2.5° and 0°, 
respectively. There is a clear gradation of properties 
and evolutionary patterns along the four columns of 
the figure. The amount of flux remaining in the head 
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Fig. 15. — Percentage of the initial magnetic flux re- 
maining in the head of the tube once the terminal ve- 
locity has been reached (stars) as a function of '^mt- 
Shown is also the percentage of flux that has been 
dragged to the wake (circles). 

of the tube is an increasing function of '^mt- This can 
be seen in Fig. 15: the stars, representing the per- 
centage of the initial flux which remain in the head, 
go from more than 80% for "ifmt = 13.9° to about 
20% for ^rnt = 2.5° and (trivially) 0% for = 0°. 
In this figure, the change between the two regimes 
occurs between, say, 3° and 8°. 

Not only the size, but also the shape of the re- 
maining head varies markedly from one case to the 
other (Fig. 14, upper panels). For \l/mt = 13.9°, the 
tube maintains its roundish shape all along. This is 
specially striking at the rear of the tube: the rising 
return flow between the two rolls in the wake cannot 
dent the back of the tube. This case is closest to the 
motion of a rigid cylinder in a fluid. For ^„it decreas- 
ing down to 2.6°, the back of the tube becomes first 
fiat and then concave through the action of the return 
flow. This transition coincides with the increase of the 
magnetic Weber number up across unity: We is 0.19, 
0.74 and 5.75 for = 13.9°, 7° and 2.5°, respec- 
tively. In the untwisted case, there is no longer a neat 
separation between tube interior and external flows. 
The rolls in the wake now contain 70% of the original 
magnetic flux, with most of the rest contained in the 
upper arch linking them. The evolutionary pattern in 
this case is similar in many respects to the results of 
Longcope et al. (1996). 

Further insight into the structural differences for 
the tubes of Fig. 14 can be obtained by comparing 
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Fig. 16. — Initial and asymptotic profiles of the dis- 
tribution of magnetic energy in the transverse field, 
^mag^ ' ^long the vertical axis of symmetry compared, 
in both cases, with the asymptotic kinetic energy den- 
sity at the apex of the tube. The asymptotic val- 
ues were calculated att = 5. The different curves cor- 
respond to cases with initial pitch angle '^mt = 13.9° 
(dashed), "^rnt = 7° (sohd) and "^rnt = 2.5° (dotted). 



the energy density of the transverse field, e^^g with 
the kinetic energy reached asymptotically by the head 
of the tube. The order of magnitude of the latter is 
(B/^/Stt) ^ (1 + tan^ *). This varies fittle (about 6% 
only) between = 2.5° and "^rnt = 13.9°. The 
magnetic energy of the transverse field, in contrast, 
varies by a large factor (approx. 30) between those 
two extremes. For the comparison (Fig. 16), we plot 



now the distribution of e^„„ along the vertical axis 
at times t = Q (upper panel) and t = 5 (lower panel), 
divided in both cases by the kinetic energy density at 
the upper stagnation point of the tube in the asymp- 
totic regime, e^^^ (calculated in the figure at time 
f = 5). In the case with the highest ^mt (dashed 
line), at t = only a thin skin at the outermost tube 
boundary is below the terminal kinetic energy level; 
we expect a skin of roughly that size to be dragged to 
the wake. For ^rnt = 7° (solid Hue) and = 2.5° 
(dotted line), most (or all) of the tube at time t = Q 
has e^„„ below the ej, level. However, in the initial 
phases of the evolution the transverse field is substan- 
tially intensified through the compression and stretch- 
ing phenomena explained in §4. This allows the head 
of the tube to be formed and its energy e^^^^ to be 
at (or above) the e^^.^^ level in the asymptotic regime 
(Fig. 16, dotted line in the second plot). 

8. Discussion 

The results presented in this paper can be dis- 
cussed from a twofold perspective. First, they can 
be applied to the theory of the magnetic activity in 
the Sun, trying to understand different aspects of the 
rise of magnetic flux across the convection zone. On 
the other hand, the simulations described here bear 
a strong resemblance to the results of laboratory and 
numerical experiments on the motion of air bubbles 
and rigid cylinders; it would be as well to clarify simi- 
larities and differences. The following subsections dis- 
cuss different aspects of those topics. 

8.1. The parameters and distribution of bcisic 

physical variables 

The numerical calculation of the present paper is 
necessarily idealized. The calculations were done in 
a high-beta regime with small ratio of radius to pres- 
sure scaleheight. In thus far, they are within the pa- 
rameter regime of expected for the magnetic tubes in 
the deep convection zone. The latter have a larger /3, 
by a factor between 10 and 100, than adopted here 
but many qualitative features of the rise should be 
similar in both cases. The only diffusive process con- 
sidered, though, the ohmic resistivity, together with 
the numerical viscosity and diffusion, yield laminar 
flows in and around the tube with Reynolds number 
of a few hundred. This cannot be expected to hold 
in the actual Sun: there, the rise of the tube will be 
accompanied by turbulent flows (with, in particular. 



21 



a turbulent wake). 

In spite of the foregoing, a conclusion seems un- 
avoidable: for the transport of the magnetic flux to 
the surface in the form of buoyant magnetic flux tubes 
to be effective, the latter must be twisted from the 
early stages of their rise. We have shown for a tube 
with a density deficit of order the full isothermal value 
{Ap/p ~ — 1//3) that an average pitch angle around 
sin~^[(i?/7fp)^/^] is indeed necessary for the tube to 
withstand the various deforming agents. Otherwise, 
they lose most of their magnetic flux to the trailing 
wake. As shown by earlier authors, the vortex rolls of 
the wake can easily end up moving horizontally rather 
than rising. The condition on the minimum pitch an- 
gle may be less stringent, though, if we consider that 
the rise may actually be driven by the Parker insta- 
bility in tubes which are stored in a neutrally buoyant 
equilibrium (Caligari et al. 1995). The instability oc- 
curs with the upgoing mass elements being driven by 
a vertical force which remains well below pg/P while 
the amplitude is not too large. To a driving force 
of amount f pg/P (/ <C 1) corresponds through (11) 
a minimum pitch angle which is about a factor /^/^ 
smaller than for isothermal straight tubes. For exam- 
ple, in the case discussed in the paper of Caligari et 
al. 1995, / < 0.04 for the top of the rising loop in the 
lowermost 10,000 km of the rise, i.e., for distances 
equivalent to about 5 times the tube radius. 

Even if the tubes are strongly buoyant [|Ap//7| = 
0(l//3) ] from the beginning, their time evolution may 
depend to some extent on the precise distribution of 
the density deficit in the tube at time t = 0. A gaus- 
sian profile, as used in this paper, yields a fast defor- 
mation of the initial magnetic configuration, followed 
by internal torsional oscillations. This is directly asso- 
ciated with the torque of the buoyancy force through- 
out the tube interior. A top-hat profile (pi < pg 
but Ap =const), on the other hand, has no buoyancy 
torques associated except at the tube boundary. Yet, 
the time evolution of such a tube for smooth pitch an- 
gle distributions like (8) or (9) is qualitatively similar 
to the cases studied in this paper, as has been tested 
through a series of numerical experiments. This can 
be understood in different ways: first, the hydrody- 
namic forces associated with the external flow depend 
on the speed of rise of the tube, which is a function of 
the average buoyancy rather than of the precise shape 
of the Ap profile. To counteract them, one needs to 
have a sufficient pitch angle in the tube interior. On 
the other hand, the gravitational term in the vortic- 



ity equation (14) for a top-hat profile is c;oncentrated 
at the periphery of the Ap distribution, but is corre- 
spondingly more intense. To counteract it, it is nec- 
essary to have a large current density J^, which can 
only be achieved if Bf is suflaciently intense, at least 
close to the tube boundary. Through both arguments 
one arrives at criteria basically like Eq. (11). 

8.2. Twisted tubes in the convection zone: 
three-dimensional effects 

If the magnetic tubes that yield active regions at 
the surface must be twisted already in the early stages 
of the rise, the Sun should have a mechanism to rou- 
tinely produce twisted tubes in the dynamo layers. 
If one relaxes the condition of symmetry along the 
tube imposed in the present paper, one may think 
that the twist could result through torsional shear- 
ing. This might come about as a natural subproduct 
of the generation of vorticity explained in this pa- 
per, if it took place at different rates on the different 
cross sections of the tube along the axis. This can 
be the case, for instance, if the originally horizontal 
tube develops an omega-loop shape in which a stretch 
of it rises while the rest remains at the original level. 
If the tube was originally untwisted (or only weakly 
twisted), the rising sections would tend to turn into 
vortex tube pairs which contain most of the original 
magnetic flux. Now, in contrast to the 2D case, this 
rotation produces a transverse fleld component in the 
flanks of the rising section of the loop. If enough Bt is 
built up, the rotation in the vortices could be braked. 
Yet, a simple calculation shows that for this mecha- 
nism to be at all effective, the footpoint separation, 
A, of the omega loop has to be small, in fact smaller 
than one pressure scaleheight. More precisely, if we 
require that the tube develops a transverse compo- 
nent of the level required by the criterion (11) after 
rising a distance equivalent to a few of its own radii, 
we obtain 

This is a very small footpoint separation: the ratio 
of the Lorentz force associated with the curvature of 
the field lines to the buoyancy force in such a tube is 
at least of order Hp/\. Consequently, it is not easy 
to raise a loop that is narrow enough for a sufficient 
pitch angle to develop at all. This is an indication 
that the tubes possibly have to be formed with the 
necessary amount of twist before they start to rise. 
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An interesting possibility for generating magnetic 
tubes with a non-zero total twist has been discussed 
by Cattaneo, Chiueh, & Hughes 1990: they consider 
a Ray leigh- Taylor unstable slab of plasma with hor- 
izontal magnetic field such that the angle subtended 
by the field vector with a fixed horizontal axis mono- 
tonically changes with depth. Upon development of 
the instability, the layer yields magnetic tubes with 
non-zero degree of twist. In their experiment, the au- 
thors showed how the resulting tubes were more re- 
sistant to deformation through the surrounding flows 
than in the corresponding untwisted case (Cattaneo 
& Hughes 1988). A corresponding three-dimensional 
problem has been calculated by Matthews, Hughes & 
Proctor 1995. They show how an unstable layer with 
a parallel horizontal field produces magnetic tubes 
with non-zero vorticity. Through nonlinear interac- 
tion, these tubes arch as they rise in a vertical plane, 
thereby becoming twisted. 

Another open question concerns the fate of a twisted 
tube as it rises across successive density scaleheights 
in the convection zone. The internal and external 
densities are essentially equal for most of the rise (ex- 
cept for the small relative difference of order up to 
I IP 1): the tube thus expands and its field in- 
tensity weakens by several orders of magnitude along 
the journey (see Moreno-Insertis 1986, 1992, 1997b). 
In fact, the beginning of this process can be seen by 
comparing the max;imum field intensity of the three 
panels of Fig. 8. The rate of decrease associated with 
an off-axis expansion of the tube, though, is different 
for the longitudinal and transverse components of the 
field (Parker 1979, §9). A rough but simple argument 
based on the conservation of fiux in 2D yields a rate 
of change of the twist following an approximate law 
tan ^ (x R. If so, the tubes could reach the upper con- 
vection zone with very large degrees of twist, possibly 
such that they can become kink unstable. Yet, three- 
dimensional effects could render that simple law of 
little use. The stretching of the tube apex in the lon- 
gitudinal direction, for instance, a phenomenon com- 
mon in Parker-unstable rising loops, may reduce the 
level of twist there. Also, the conservation of the to- 
tal helicity in the tube could put an upper limit to 
the number of turns of the field lines around the axis 
at any single place. However it is, large pitch an- 
gles should be a common appearance in rising tubes 
when they reach photospheric levels. The sheared 
field structures observed in emerging active regions 
(Lites et al.l995, Leka et al.l996) may be a conse- 



quence of this (sec also Tanaka 1991, Kurokawa 1989 
and Rust & Kumar 1996). 

8.3. Magnetic tubes, air bubbles and rigid 
cylinders: similarities and differences 

Along this paper we have compared our results to 
the laboratory experiments of flow past rigid cylinders 
and air bubbles at Re = 200 300. In the following 
we summarize the similarities and differences found. 

Magnetic tubes and air bubbles. Remarkable simi- 
larities are: (a) the general shape and structural fea- 
tures (skirt, central tail) (b) the protection of the 
interior of the rising object through surface tension. 
These similarities are all the more striking given the 
difference in buoyancy [density deficit \Ap/p\ <C 1 for 
the tubes, 0(1) for the bubbles] and in the physi- 
cal source of the surface tension (capillary effects vs. 
jump in the tangential field component). As a result, 
defining the Weber number on the basis of the corre- 
sponding surface tension mechanism, one can formu- 
late a common law of dependence of structural prop- 
erties on We. 

Magnetic tubes and rigid cylinders. A clear similar- 
ity concerns the boundary layer: both objects have a 
no-slip condition along the boundary. Consequently, 
the relative jump of the tangential velocity across the 
boundary layer is large. [Vig]/v^g ^ 0(1)- In con- 
trast, there is a zero-tangential stress condition along 
the fluid/ air interface of a bubble; correspondingly, 
[vtgl/Vfg ~ 0(i?e~^/^). As a result of the no-sHp con- 
dition, both buoyant magnetic tubes and rigid cylin- 
ders produce (and shed downstream) secondary rolls 
near the point of separation of the boundary layer. 

Features specific to the buoyant magnetic tubes with 
twist. In the magnetic tubes of the present simu- 
lations, the wake is formed in the initial phases of 
the rise out of material bodily transported from the 
initial tube Thus, the wake is magnetized and main- 
tains a magnetic connection to the head of the rising 
tube along time. The kind of behavior followed by 
the magnetic rope depends directly on the amount of 
magnetic flux incorporated into the wake, which, in 
turn, depends on the initial twist. In the case of an 
initially highly twisted tube {^mt = 13.9°), the mag- 
netic rope is almost rigid and its shape and the flow 
around it strongly remind those of a solid cylinder. If 
the initial pitch angle is closer to the threshold (11), 
then the tube deforms and adopt a bubble-like shape 
(though it still satisfies a no-slip condition along its 
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boundary). Finally, if "^mt is very small, the tube 
behaves like a rising thermal (Longcope et al. 1996). 
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